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We show that massive lepton-pair production, the Drell-Yan process, should be 
a good source of independent constraints on the polarized gluon density, free from 
the experimental and theoretical complications of photon isolation that beset stud- 
ies of prompt photon production. We provide predictions for the spin-averaged and 
spin-dependent differential cross sections as a function of transverse momentum at 
energies relevant for the Relativistic Heavy Ion Collider (RHIC) at Brookhaven. 
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1. Introduction 



Massive lepton-pair production, hi + h 2 — > 7* + X; 7* — > /f, known as the Drell-Yan 
process *\ and prompt real photon production, h± + h 2 — > 7 + X, are two of the most 
valuable probes of short- distance behavior in hadron reactions. They supply critical infor- 
mation on parton momentum densities and opportunities for tests of perturbative quan- 
tum chromodynamics (QCD). Spin- averaged parton momentum densities may be extracted 
from spin-averaged nucleon-nucleon reactions, and spin-dependent parton momentum den- 
sities from spin-dependent nucleon-nucleon reactions. An ambitious experimental program 
of measurements of spin-dependence in polarized proton-proton reactions will begin soon 
at Brookhaven's Relativistic Heavy Ion Collider (RHIC) with kinematic coverage extending 
well into the regions of phase space in which perturbative quantum chromodynamics should 
yield reliable predictions. 

The Drell-Yan process has tended to be thought of primarily as a source of information on 
quark densities. Indeed, the mass and longitudinal momentum (or rapidity) dependences of 
the cross section (integrated over the transverse momentum Qt of the pair) provide essential 
constraints on the antiquark momentum density, complementary to deep-inelastic lepton 
scattering from which one gains information of the sum of the quark and antiquark densities. 
Prompt real photon production, on the other hand, is a source of essential information on the 
gluon momentum density. At lowest order in perturbation theory, the reaction is dominated 
at large values of the transverse momentum pr of the produced photon by the "Compton" 
subprocess, q + g — > 7 + q. This dominance is preserved at higher orders, indicating that the 
experimental inclusive cross section differential in p T may be used to determine the density 
of gluons in the initial hadrons. 2 ) 

In this paper, we summarize recent work 3 )' 4 ), in which we demonstrate that the Comp- 
ton subprocess, q + g — > 7* + q also dominates the Drell-Yan cross section in polarized 
and unpolarized proton-proton reactions for values of the transverse momentum Qt of the 
pair that are larger than roughly half of the pair mass Q, Qt > Q/2. The Drell-Yan pro- 
cess is therefore a valuable, heretofore overlooked, independent source of constraints on the 
spin-averaged and spin-dependent gluon densities. Although the Drell-Yan cross section is 
smaller than the prompt photon cross section, massive lepton pair production is cleaner the- 
oretically since long-range fragmentation contributions are absent as are the experimental 
and theoretical complications associated with isolation of the real photon. 5 ) Moreover, the 
dynamics of spin-dependence in hard-scattering processes is a sufficiently complex topic, and 
its understanding at an early stage in its development, that several defensible approaches 
for extracting polarized parton densities deserve to be pursued with the expectation that 
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consistent results must emerge. 
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In this Section differential cross sections for 
massive lepton-pair production are presented as 
functions of Qt at collider energies. We work in the 
MS renormalization scheme and set the renormal- 
ization and factorization scales equal. We employ 
the MRST98-1 set of spin-averaged parton densi- 
ties 6 -* and a two-loop expression for the strong cou- 
pling strength a s (fj), with five flavors and appro- 
priate threshold behavior at /x = m b . We use the 
value = 300 MeV of the MRST98-1 set. The 
strong coupling strength a s is evaluated at a hard 
scale n = \JQ 2 + Qt- We offer results for three 
values of RHIC center-of-mass energy, v^S* = 50, 
200, and 500 GeV. 

For v^S* = 200 GeV, we present the invariant 
inclusive cross section Ed 3 a /dp 3 as a function of 
Qt in Fig. 1. Shown are the qq and qg perturba- 
tive contributions to the cross section at leading 
order and at next-to- leading order. For Qt < 1.5 
GeV, the qq contribution exceeds that of qg chan- 
nel. However, for values of Q T > 1.5 GeV, the qg 
contribution becomes increasingly important. The 
qg contribution accounts for about 80 % of the rate 
once Qt — Q. Subprocesses other than those ini- 
tiated by the qq and qg initial channels contribute 
negligibly. 

In Ref. 4), results are shown also for a larger value of Q: 11 < Q < 12 GeV, and for the 
interval 2.0 < Q < 3.0 GeV. The fractions of the cross sections attributable to qg initiated 
subprocesses again increase with Q Tl growing to 80 % for Q T ^ Q. In region 2.0 < Q < 3.0 
GeV, one would doubt the reliability of leading-twist perturbative descriptions of the cross 
section da/dQ, integrated over all Qt- However, for values of Qt that are large enough, 
cx s {Qt) is small and a perturbative description of the Qt dependence of d?a / dQdQT ought 



Fig. 1. Lowest order and next-to- 
leading order perturbative calcula- 
tions of Ed 3 a /dp 3 as a function 
of Q T for pp — > 7*V at \fS = 
200 GeV. Contributions from the 
qg and qq channels are shown sep- 
arately. The results are averaged 
over the rapidity interval - 1.0 < 
y < 1.0 and over the interval 5.0 
< Q < 6.0 GeV. 
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to be justified. 

In Fig. 2, we provide next-to-leading order pre- 
dictions of the differential cross section func- 
tion of Qt for three values of the center-of-mass 
energy. In order to ascertain the range of values 
of Qt that can be explored, we take Ed 3 a/dp 3 = 
l(T 3 pb/GeV 2 as the minimum accessible cross sec- 
tion. This level is based on assumed luminosi- 
ties 7 ) of 8 x lO^cm^sec -1 at VS = 200 GeV, and 
2 x lO^cm^sec^ 1 at = 500 GeV, along with 
runs of 10 weeks per year, equivalent to integrated 
luminosities of 320 pb _1 year _1 and 800 pb _1 year _1 . 
Better performance than anticipated of the accel- 
erator and extended running time will increase the 
reach. The luminosity scales roughly with v^S* be- 
cause the focusing power can be increased as the 
energy is reduced. This capability saturates near 
70 GeV/c per beam, and below about y/S = 140 
GeV the luminosity drops roughly as S. 

Adopting the nominal value 
Ed 3 a/dp 3 = 10~ 3 pb/GeV 2 , we use the curves in 
Fig. 2 to establish that the massive lepton-pair 
cross section may be measured to Qt = 7.5, 14, and 
18.5 GeV at = 50, 200, and 500 GeV, respectively, when 2 < Q < 3 GeV, and to Q T = 6, 
11.5, and 15 GeV when 5 < Q < 6 GeV. In terms of reach in the fractional momentum x g \ uon 
carried by the gluon, these values of Qt may be converted to x g i uon ~ xt = 2Q T /VS = 0.3, 
0.14, and 0.075 at y/S = 50, 200, and 500 GeV when 2 < Q < 3 GeV, and to x gluon ~ 0.24, 
0.115, and 0.06 when 5 < Q < 6 GeV. On the face of it, the smallest value of y/S provides 
the greatest reach in x g i uon . However, the reliability of fixed-order perturbative QCD as well 
as dominance of the qg subprocess improve with greater Q T . The maximum value Qt — 7.5 
GeV attainable at = 50 GeV argues for a larger y/S. 

Comparing the magnitudes of the prompt photon and massive lepton pair production 
cross sections 4 \ we note that the inclusive prompt photon cross section is a factor of 1000 to 
4000 greater than the massive lepton-pair cross section integrated over the mass interval 2.0 
< Q < 3.0 GeV, depending on the value of Q T - This factor is attributable in large measure 
to the factor a em /(3nQ 2 ) associated with the decay of the virtual photon to Again 
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Fig. 2. Differential cross sections Ed 3 a /dp 3 
function of Qt for for pp — > 
7* + X at ^/S = 50, 200, and 500 
GeV, averaged over the rapidity in- 
terval -1.0 < y < 1.0 and the mass 
interval 5.0 < Q < 6.0 GeV. 
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taking Ed 3 a/dp 3 = 10 _3 pb/GeV as the minimum accessible cross section, we establish that 
the real photon cross section may be measured to pt = 14, 33, and 52 GeV at v^S* = 50, 
200, and 500 GeV, respectively. The corresponding reach in xt is xt = 2p T /\fS = 0.56, 
0.33, and 0.21 at VS = 50, 200, and 500 GeV. 

The significantly smaller cross section in the case of massive lepton-pair production means 
that the reach in x g i uon is restricted to a factor of about two to three less, depending on 
and Q, than that potentially accessible with prompt photons in the same sample of data. 
Nevertheless, it is valuable to be able to investigate the gluon density with a process that 
has reduced experimental and theoretical systematic uncertainties from those of the prompt 
photon case. 

In Ref. 3) we compared our spin-averaged cross sections with available fixed-target and 
collider data on massive lepton-pair production at large values of Qt, and we were able to 
establish that fixed-order perturbative calculations, without resummation, should be reliable 
for Qt > Q/2. The region of small Qt and the matching region of intermediate Qt are com- 
plicated by some level of phenomenological ambiguity. Within the resummation approach, 
phenomenological non-perturbative functions play a key role in fixing the shape of the Qt 
spectrum at very small Qt, and matching methods in the intermediate region are hardly 
unique. For the goals we have in mind, it would appear best to restrict attention to the 
region Q T > Q/2. 



§3. Predictions for Spin Dependence 

Given theoretical expressions that relate the spin-dependent cross section at the hadron 
level to spin-dependent partonic hard-scattering matrix elements and polarized parton densi- 
ties, we must adopt models for spin-dependent parton densities in order to obtain illustrative 
numerical expectations. The current deep inelastic scattering data do not constrain the po- 
larized gluon density tightly, and most groups present more than one plausible parametriza- 
tion. Gehrmann and Stirling (GS) 8 ) present three such parametrizations, labelled GSA, 
GSB, and GSC. In the GSA and GSB sets, AG(x,/j, ) is positive for all x, whereas in the 
GSC set AG(x,/jl ) changes sign. After evolution to fij = 100 GeV 2 , AG(x,fif) remains 
positive for essentially all x in all three sets, but its magnitude is small in the GSB and 
GSC sets. We use the three parametrizations suggested by Gehrmann and Stirling, and we 
verified that the positivity requirement Afl(x, /if)/ fl(x, /if) < 1 is satisfied. 

Other sets of spin-dependent parton densities have been published, e.g., the set of Gliick, 
Reya, Stratmann, and Vogelsang (GRSV). 9 ) The three GS parametrizations of the polarized 
gluon density span a range of possibilities that is very similar to that spanned by the four 
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gluon densities of GRSV. It is not our purpose to "test", or to suggest tests, of existing 
parametrizations all of which will have been modified substantially by the time data are 
available from RHIC. Rather, we use the existing parametrizations as illustrative possibilities 
in order to estimate the range of magnitudes that might be expected for All at RHIC energies 
and to gauge the sensitivity that measurements of spin dependence in massive lepton-pair 
production may offer. The goal, after all, is to measure the polarized gluon density, not to 
test parametrizations. 

In this Section, we present two-spin longitudi- 
nal asymmetries for massive lepton-pair production 
as a function of transverse momentum. Results are 
displayed for pp collisions at the center-of-mass en- 
ergies y/S = 50, 200, and 500 GeV typical of the 
Brookhaven RHIC collider. 

The two-spin longitudinal asymmetries, All, 
are computed in leading-order. More specifi- 
cally, we use leading-order spin-averaged and spin- 
dependent partonic subprocess cross sections, a 
and Aa, with next-to-leading order spin-averaged 
and spin-dependent parton densities and a two- 
loop expression for a s . The choice of a leading- 
order expression for Aa is required because the full 
next-to- leading order derivation of Aa has not been 
completed for massive lepton-pair production. Ex- 
perience with prompt photon production indicates 
that the leading-order and next-to-leading order re- 
sults for the asymmetry are similar so long as both 
are dominated by the qg subprocess. To obtain the 
spin-dependent cross section, we use the GS polar- 
ized parton densities with the GS value = 231 
MeV. This value of A^ differs from the MRST 
value A^ = 300 MeV for the unpolarized densities 
used in our computation of the spin-averaged cross 

section. The use of different values of A^ for the spin-dependent and spin-averaged cross 
sections may appear unfortunate. However, there is not much of an alternative at present 
short of creating new sets of polarized parton densities based on the most up-to-date spin- 
averaged densities that we prefer to use for the spin-averaged cross section. To change A^ 



Qt (GeV) 

Fig. 3. Computed longitudinal asym- 
metry All as a function of Qt for 
pp -> ~/*X at = 200 GeV. The 
asymmetry is averaged over the ra- 
pidity interval -1.0 < y < 1.0 and 
mass interval 5.0 < Q < 6.0 GeV, 
for three choices of the polarized 
parton densities, GSA, GSB, and 
GSC. 
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arbitrarily in the GS set to equal that of the MRST set would distort the spin-dependent 
densities. 

In Fig. 3, we present the two-spin longitudinal asymmetries, All, as a function of Qt for 
the three GS choices of the polarized gluon density. The asymmetry becomes sizable for large 
enough Qt for the GSA and GSB parton sets but not in the GSC case. The asymmetry All 
is nearly independent of the pair mass Q as long as Qt is not too small. This feature should 
be helpful for the accumulation of statistics; small bin-widths in mass are not necessary, but 
the J/ip and T resonance regions should be excluded. 

As noted above the qg subprocess dominates the 
spin-averaged cross section. It is interesting and 
important to inquire whether this dominance per- 
sists in the spin-dependent situation. In Fig. 4, we 
compare the contribution to the asymmetry from 
the polarized qg subprocess with the complete an- 
swer for all three sets of parton densities. The qg 
contribution is more positive than the full answer 
for values of Qt that are not too small; the full an- 
swer is reduced by the negative contribution from 
the qq subprocess for which the parton-level asym- 
metry oll = — 1- At small Qt, the net asymme- 
try may be driven negative by the qq contribution, 
and based on our experience with other calcula- 
tions, from processes such as gg that contribute in 
next-to-leading order. For the GSA and GSB sets, 
we see that once it becomes sizable (e.g., 5% or 
more), the total asymmetry from all subprocesses 
is dominated by the large contribution from the qg 
subprocess. 

As a general rule in studies of polarization phe- 
nomena, many subprocesses can contribute small 
and conflicting asymmetries. Asymmetries are 
readily interpretable only in situations where the 
basic dynamics is dominated by one major subpro- 
cess and the overall asymmetry is sufficiently large. 
In the case of massive lepton-pair production that 

is the topic of this paper, when the overall asymmetry All itself is small, the contribution 
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Fig. 4. Comparison of the contribution 
of the qg subprocess (upper curve) 
to the longitudinal asymmetry All 
with the total (lower curve) as a 
function of Qt for pp — > j*X at 
\/~S = 200 GeV. The asymmetry is 
averaged over the rapidity interval 
-1.0 < y < 1.0 and over the inter- 
val 5.0 < Q < 6.0 GeV. Results 
are shown for three sets of spin- 
dependent parton densities. 
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from the qg subprocess cannot be said to dominate the answer. However, if a large asym- 
metry is measured, similar to that expected in the GSA case at the larger values of Qt, 
Fig. 4 shows that the answer is dominated by the qg contribution, and data will serve to 
constrain AG(x, fXf). If AG(x,/if) is small and a small asymmetry is measured, such as for 
the GSC parton set, or at small Qt for all parton sets, one will not be able to conclude 
which of the subprocesses is principally responsible, and no information could be adduced 
about AG(x, ///), except that it is small. 

In Fig. 5, we examine energy dependence. For 
Qt not too small, All in massive lepton pair 
production is well described by a scaling func- 
tion of x T = 2Q T /y/S, A LL (y/S,Q T ) ^ h r (x T ). 
In our discussion of the spin-averaged cross sec- 
tions, we took Ed 3 a/dp 3 = 10" 3 pb/GeV as the 
minimum accessible cross section. Longitudinal 
asymmetries All = 20%, 7.5%, and 3% are pre- 
dicted at this level of cross section at VS = 50, 
200, and 500 GeV when 2 < Q < 3 GeV, and 
A LL = 11%, 5%, and 2% when 5 < Q < 6 GeV. 
For a given value of Qt, smaller values of \fS re- 
sult in greater asymmetries because AG(x)/G(x) 
grows with x. 

The predicted cross sections in Fig. 2 and the 
predicted asymmetries in Fig. 5 should make it pos- 
sible to optimize the choice of center-of-mass en- 
ergy at which measurements might be carried out. 
At v^S* = 500 GeV, asymmetries are not apprecia- 
ble in the interval of Qt in which event rates are 
appreciable. At the other extreme, the choice of 
= 50 GeV does not allow a sufficient range 
in Qt- Accelerator physics considerations favor 
higher energies since the instantaneous luminosity 
increases with a/S*. Investigations in the energy in- 
terval = 150 to 200 GeV would seem preferred. 

So long as Q T > Q, the asymmetry in massive 
lepton-pair production is about the same size as that in prompt real photon production, as 
might be expected from the strong similarity of the production dynamics in the two cases. 



upper: only qg channel 
lower: all contributions 
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Fig. 5. Computed longitudinal asym- 
metry All as a function of Qt for 
for pp -► 7* + X at y/S = 50, 
200, and 500 GeV, averaged over 
the rapidity interval -1.0 < y < 1.0 
and the mass interval 5.0 < Q < 
6.0 GeV. Shown are both the com- 
plete answer at leading-order and 
the contribution from the qg sub- 
process. The GSA set of polarized 
parton densities is used. 
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As in massive lepton-pair production, All in prompt photon production is well described by 
a scaling function of xt = 2pt/VS, All(v / 5',Pt) — h^(xx)- For Ed 3 a jdp 3 = lCT 3 pb/GeV 2 , 
we predict longitudinal asymmetries All = 31%, 17%, and 10% in real prompt photon 
production at \/~S = 50, 200, and 500 GeV. 

Although the qg Compton subprocess is dominant, one might question whether uncer- 
tainties associated with the quark density compromise the possibility to determine the gluon 
density. In this context, it is useful to recall 10 ^ that when the Compton subprocess is domi- 
nant, the spin- averaged cross section may be rewritten in a form in which the quark densities 
do not appear explicitly: 



- ~ J d Xl Jdx 2 — G(x2,Vf)—j-p— + ( Xl ~ x 2 ) . (3-1, 



dp 3 J J \ xi 1 dp 3 

Likewise, the spin-dependent cross section may be expressed as 10 ) 

Ed3 ^ lh2 « / dxi J dx 2 (2gi(xi, fi)AG(x 2 , fi)^^ + (*i <- ^ ■ (3-2) 

In Eqs. (3-1) and (3-2), F 2 (x, /Aj) and gi(x, fj/j) are the proton structure functions measured in 
spin-averaged and spin-dependent deep-inelastic lepton-proton scattering. It is evident that 
the production of massive lepton-pairs at large enough Qt will determine the gluon density 
provided the proton structure functions are measured well in deep-inelastic lepton-proton 
scattering. 



§4. Discussion and Conclusions 

In this paper we summarize a calculation of the longitudinal spin-dependence of massive 
lepton-pair production at large values of transverse momentum. We provide polarization 
asymmetries as functions of transverse momenta that may be useful for estimating the fea- 
sibility of measurements of spin-dependent cross sections in future experiments at RHIC 
collider energies. The Compton subprocess dominates the dynamics in longitudinally po- 
larized proton-proton reactions as long as the polarized gluon density AG(x,/if) is not too 
small. As a result, two-spin measurements in massive lepton-pair production in polarized pp 
scattering should constrain the size, sign, and Bjorken x dependence of AG(x,/if). 

Significant values of All (i-e., greater than 5 %) may be expected for xt = 2Qt/V~S > 
0.10 if the polarized gluon density AG(x,/j,f) is as large as that in the GSA set of polarized 
parton densities. If so, the data could be used to determine the polarization of the gluon 
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density in the nucleon. On the other hand, for small AG(x, ///), dominance of the qg 
subprocess is lost, and AG(x,fif) is inaccessible. 

Various methods have been discussed in the literature to access the spin-averaged and 
spin-dependent gluon densities. These include inclusive/isolated prompt photon production 
at large transverse momentum (with or without a tagged recoil jet); the Drell-Yan process at 
large Qt, as advocated in this paper; hadronic jet production at large pt] and heavy flavor 
production such as charm (c) and bottom (b), thought to be mediated by the subprocess 
g + g — > c + c + X. Each has its possibilities and drawbacks. Hadronic jet production 
benefits from the largest rate, but the large number of contributing subprocesses and the 
complications of jet definition mitigate against it. Inclusive prompt photon production 
is theoretically clean but for the non-perturbative long-range fragmentation contribution. 
However, experimenters measure isolated photons, and isolation renders the contact with 
theory somewhat murky. 5 ) The Drell-Yan process seems ideal, but its rate is relatively low. 
Heavy flavor production is an enigma. The charm quark is so light that calculations based 
on fixed-order perturbative QCD are of questionable reliability at collider energies. The 
bottom quark seems to be heavy enough to justify perturbation theory, but why does the 
experimental cross section at Fermilab collider energies exceed next-to-leading order QCD 
predictions by a factor of 2 or 3? 

The dynamics of spin-dependence in hard-scattering processes is a sufficiently complex 
topic, and its understanding at an early stage in its development, that several defensible 
approaches for extracting polarized parton densities deserve to be pursued with the expec- 
tation that consistent results must emerge. For the first time, RHIC offers an opportunity 
to explore the dynamics of spin dependence in hadron reactions at values of the hard scale 
that are sufficiently large that perturbative QCD ought to work. It will take several years of 
experimentation and analysis for a clear and consistent picture to emerge of spin-dependent 
parton densities, just as it has in the study of spin- averaged parton densities. The long term 
support of this program by the Brookhaven Laboratory and by the funding agencies will be 
essential for realization of the considerable potential that RHIC offers. 

Acknowledgments 

ELB is pleased to acknowledge the warm hospitality of the RIKEN Laboratory and of 
the Symposium hosts, particularly Shunzo Kumano, Masayasu Ishihara, Shunichi Kobayashi, 
Toshi-Aki Shibata, and Koichi Yazaki. Work in the High Energy Physics Division at Ar- 
gonne National Laboratory is supported by the U.S. Department of Energy, Division of 
High Energy Physics, Contract W-31-109-ENG-38. This work was supported in part by 



10 



DOE contract DE-AC05-84ER40150 under which Southeastern Universities Research Asso- 
ciation operates the Thomas Jefferson National Accelerator Facility. MK is supported by 
Bundesministerium fur Bildung und Forschung under Contract 05 HT9GUA 3, by Deutsche 
Forschungsgemeinschaft under Contract KL 1266/1-1, and by the European Commission 
under Contract ERBFMRXCT980194. 

References 

1) S. Drell and T. M. Yan, Phys. Rev. Lett. 25 (1970), 316; Ann. Phys. (NY) 66 (1971), 
578. 

2) E. L. Berger and J.-W. Qiu, Phys. Lett. B248 (1990), 371 and Phys. Rev. D44 
(1991), 2002; H. Baer, J. Ohnemus, and J. F. Owens, Phys. Rev. D42 (1990), 61; P. 
Aurenche et al, Nucl. Phys. B399 (1993), 34; L. E. Gordon and W. Vogelsang, Phys. 
Rev. D48 (1993), 3136 and D50 (1994), 1901; M. Gliick, L. E. Gordon, E. Reya, 
and W. Vogelsang, Phys. Rev. Lett. 73 (1994), 388; L. E. Gordon, Nucl. Phys. B501 

(1997) , 175. 

3) E. L. Berger, L. E. Gordon, and M. Klasen, Phys. Rev. D58 (1998), 074012 and 
hep-ph/9906402. 

4) E. L. Berger, L. E. Gordon, and M. Klasen, hep-ph/9909446, to be published in 
Phys. Rev. D. 

5) E. L. Berger, X.-F. Guo, and J.-W. Qiu, Phys. Rev. Lett. 76 (1996), 2234; Phys. 
Rev. D54 (1996), 5470; and in '97 QCD and High Energy Hadronic Interactions, 
Proceedings of the XXXIInd Rencontres de Moriond, Les Arcs, France, March, 1997, 
edited by J. Tran Thanh Van (Editions Frontieres, Paris, 1997) pp 267 - 274. 

6) A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thorne, Eur. Phys. J. C4 

(1998) , 463; A. D. Martin, R. G. Roberts, W. J. Stirling, and R. S. Thorne, hep- 
ph/9907231. 

7) G. Bunce, Brookhaven National Laboratory, private communication. 

8) T. Gehrmann and W. J. Stirling, Phys. Rev. D53 (1996), 6100. 

9) M. Gliick, E. Reya, M. Stratmann, and W. Vogelsang, Phys. Rev. D53 (1996), 4775. 
10) E. L. Berger and J.-W. Qiu, Phys. Rev. D40 (1989), 778 and D40 (1989), 3128. 



11 



